O 
O 

(N 



X 



Geometrical Spinoptics 
and 

the Optical Hall Effect 

c. DUVAL* 

Centre de Physique Theorique, CNRS, Luminy, Case 907 
F- 13288 Marseille Cedex 9 (France) § 

Z. HORVATH 1 
Institute for Theoretical Physics, Eotvos University 
Pazmany P. setany 1/A 



> 

O : H-1117 BUDAPEST (Hungary) 

^ ; P. A. HORVATHY 11 



Laboratoire de Mathematiques et de Physique Theorique 



O 

p ,. Universite de Tours, Pare de Grandmont 

Jl ' F-37 200 TOURS (France) 

> ■ 



Abstract 

Geometrical optics is extended so as to provide a model for spinning light 
rays via the coadjoint orbits of the Euclidean group characterized by color and 
spin. This leads to a theory of "geometrical spinoptics" in refractive media. 
Symplectic scattering yields generalized Snell-Descartes laws that include the 
recently discovered optical Hall effect. 



MSC : 78A05; 37J15; 70H40 

Preprint : CPT-2005/P.052 arXiv:math-ph/0509031 

■^mailto: duval@cpt.univ-mrs.fr 

§ UMR 6207 du CNRS associee aux Universites d'Aix-Marseille I et II et Universite du Sud 
Toulon- Var; Laboratoire affilie a la FRUMAM-FR2291 
^mailto: zalanh@ludens.elte.hu 
II mailto: horvathy@lmpt.univ-tours.fr 



1 Introduction 



The mathematical foundation of optics certainly goes back to Euclid's Optics with 
the observation that light travels along straight lines, namely along the geodesies of 
Euclidean space. Geometrical optics has, since the breakthrough provided by the 
Fermat principle of least optical path, proved extremely useful as a theory of light, 
prior to Maxwellian wave optics. Suffice it to mention its relevance in the design of 
optical (as well as electronic!) devices using lenses, mirrors, etc., its importance in 
the understanding of caustics [T], in the modelling of optical aberrations, etc. Geo- 
metrical optics has, since then, been recognized as a "semi-classical" limit of wave 
optics with small parameter X (where A is a typical wavelength); it has nevertheless 
been constantly considered as a self-consistent theory for light rays, borrowing much 
from differential geometry, and, more specifically, from Riemannian and symplec- 
tic geometries. Geometrical optics provides, indeed, a beautiful link between both 
previously mentioned geometries: (i) light travels along (oriented) geodesies of an 
optical medium, a 3-dimensional manifold whose Riemannian structure is defined 
by the refractive index, (ii) the set of all such geodesies is naturally endowed with a 
structure of 4-dimensional symplectic manifold. It is this duality that will serve as 
an Ariadne's thread in our subsequent extension of geometrical optics. 

In addition, it is nowadays a well-established experimental fact that trajecto- 
ries of light beams in inhomogeneous media slightly depart from those enacted by 
the Fermat principle. This class of effects predicted by Fedorov, fifty years ago, 
has since then received numerous theoretical interpretations that go back to work 
of Costa de Beauregard Boulware among many a fundamental and more 
recent contribution (see, e.g., jH[7j for a brief historical account with an updated list 
of references regarding both theory and experiment). One the measured effects is 
the "Magnus Effect for light" that describes how trajectories depend upon the polar- 
ization state of the beam in weakly inhomogeneous media. Another related effect is 
the so-called "Optical Hall Effect" (OHE) which has lately received special attention 
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and is associated with a transverse shift of the position of a photonic wave packet 
at the interface separating two media of different refractive indices. (The OHE ac- 
tually bears strong resemblance with the Hall effect governing, in two-dimensional 
conducting samples, the electronic transport, transverse to the electric field and an 
applied external magnetic field.) Such a shift, transverse to the incidence plane, has 
already been experimentally measured by Imbert [23] in the case of total reflection. 
As to the transverse shift for partial reflection and refraction, it is currently under 
highly active investigation from both a theoretical and an experimental perspective; 
see, e.g., |H [7j and [31]. Let us also mention the recently discovered phenomenon 
of magneto-transverse light diffusion in Faraday-active dielectric media [SZ], the 
"Photonic Hall Effect" (PHE), see [39J, which is clearly a spin-induced effect. 

The need for a generalization of geometrical optics which would consistently in- 
clude polarization effects hence became mandatory. Various approaches, including 
a full-fledged computation of refraction and reflection of arbitrarily polarized Gaus- 
sian electromagnetic wave packets [7j, have been put forward. Of particular interest 
are recent extensions of geometrical optics, within a Maxwellian context, using a 
certain "Berry connection" whose curvature (in momentum space) yields a modifi- 
cation of the Fermat equations of motion for polarized light beams [23 El El EJ] ; 
see also [21 El E] A quasi-classical formula for the above-mentioned transverse 
shift of polarized light beams has also been proposed by Onoda, Murakami, and 
Nagaosa [21] , together with an experimental set up using photonic crystals in order 
to reveal the OHE for reflected and refracted light beams. Nevertheless, no consen- 
sus seems to have emerged so far regarding a clear-cut theory of geometrical optics 
including polarization. 

Our standpoint is to take advantage of the fact that geometrical optics is fun- 
damentally related to Euclidean geometry since oriented straight lines (light rays) 
actually consist of specific coadjoint orbits of the Euclidean group, E(3). An ex- 
tended theory of geometrical optics should therefore be expected to emerge from 
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the same Euclidean geometry, more precisely from the consideration of the other 
E(3)-coadjoint orbits which seem to have been overlooked by physicists. The pur- 
pose of this article is therefore to view and exploit the generic E(3)-coadjoint orbits, 
carrying color and spin according to Souriau's classification as the sets of "free" 
colored and circularly polarized light rays. We then introduce a natural adaptation 
of the Fermat prescription to the spin case, inspired by the prescription of "minimal 
coupling" to a curved metric used in general relativity, see, e.g., [13 EBl EH EHj- 
This leads us in a straightforward fashion to a theory describing the trajectories of 
spinning light rays in arbitrary dielectric media. We name this theory "geometrical 
spinoptics" . 

The present article is organized as follows. 

Section |21 is devoted to a general overview of the geometry of the set of light rays 
with color and spin as coadjoint orbits of the group SE(3) of Euclidean orientation- 
preserving 1 isometries of Euclidean 3-dimensional space space. It should be high- 
lighted that the generic coadjoint orbits, with topology TS 2 , are automatically en- 
dowed with a "twisted" cotangent symplectic structure which describes spin — or 
polarization. There is, hence, no need to introduce a "Berry curvature" which is, in 
a sense, already encoded in the above-mentioned geometrical twist. The relationship 
of these Euclidean coadjoint orbits to the massless, spinning, coadjoint orbits of the 
Poincare and Galilei groups is spelled out explicitly. 

In Section 01 we generalize the Fermat prescription so as to describe the pre- 
symplectic structure of colored and spinning light rays in an arbitrary refractive 
medium characterized by a dielectric tensor defining a Riemannian metric. This 
prescription happens to be akin to the "minimal coupling" procedure used, in general 
relativity, to account for the geodesic deviation of spinning particles due to tidal 
forces generated by the interaction of spin and spacetime curvature. The general 

equations for spinoptics (j3.17J) are then derived. The special case of an isotropic 

1 Considering, here, the neutral component of the Euclidean group is a technicality which will 
be discussed and justified below. 
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and inhomogeneous medium, of refractive index n, is worked out in great detail; we 
recover, upon linearization of the foliation (|3.28J) around the value n = const., a 
system of equations of motion for polarized light rays proposed in |31| . 

Section lUgives us the opportunity to derive, in the present framework, the Snell- 
Descartes laws (J4.30)) for spinoptics, dictated by Souriau's symplectic scattering. 
We explicitly compute the form of the reflection and refraction (local) symplecto- 
morphism between "in" and "out" polarized optical states. A novel and subtle 
phenomenon naturally stems from our geometrical treatment, viz., the fact that 
scattered spinning light rays are actually shifted, transversally to the plane of in- 
cidence. The formula 1)4.33)1 we derive for this transverse shift agrees with the one 
proposed in [HI] to describe the OHE from a different standpoint. This shift is then 
analyzed for left-handed media, providing a plausible mechanical interpretation of 
the "perfectness" of superlenses |33) . 

At last, in Section we sum up the content of this article, and conclude by 
presenting several generalizations in prospect, e.g., a theory of geometrical spinoptics 
for Faraday-active media. 

A companion article provides an overview of geometrical spinoptics and 
discusses its relation to other physically oriented approaches found in the recent 
literature. 

Acknowledgment: We are indebted to K. Y. Bliokh for fruitful correspondence. 

2 Geometrical optics and the Euclidean group 

An oriented straight line, £, in Euclidean (affine) space (E 3 , (•,•)) is determined 
by its direction, i.e., a vector u e M 3 of unit length and an arbitrary point M G £. 
Should an origin, O E E 3 , be chosen, we may take the vector q = M — O to be 
orthogonal to u. The set of oriented, non parametrized, straight lines is plainly the 
smooth manifold 

M = = (q, u)el 3 xl 3 (u, u) = 1, (u, q) = 0}, (2.1) 
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i.e., the tangent bundle M. = TS 2 of the round sphere S 2 . 

Now TS 2 has been recognized [Hl| as a coadjoint orbit of the group, E(3), of 
Euclidean isometries and inherits, as such, an E(3)-invariant symplectic structure; 
see also |2lj . 

Let us briefly recall the general construction leading, in particular, to the pre- 
vious symplectic manifold. Start with the group SE(3) = SO (3) ix R 3 of orientation- 
preserving Euclidean isometries whose elements we denote g = (R, x) . Let p, = (£, p) 
be a point in se(3)* where se(3) = M 3 k IR 3 is the Lie algebra of SE(3). The coadjoint 
representation 2 of SE(3) reads then Coad s (£, p) = (R£+x.xRp, Rp) where x stands 
here for standard cross-product. Obviously C = ||p|| 2 = (p, p) and C = (£,p) are 
invariants of the coadjoint representation. If C = 0, then C" = \\£\\ 2 is an extra 
invariant. These are, in fact, the only invariants and fixing (C, C) or (C = 0, C") 
yields a single coadjoint orbit [22 | 130 } l3lj. 

Consider now, in full generality, a finite-dimensional Lie group G whose Lie 
algebra is denoted g. Fix /i G 0* and posit the following 1-form 

n7 = ^ "i? (2.2) 

on G where d stands for the (left-invariant) Maurer-Cartan 1-form of the group. 
It is a general fact that a = dw is a presymplectic 2-form on G which descends 
as the canonical Kirillov-Kostant-Souriau |2lH 12*71 Ell symplectic 2-form, u, on the 
coadjoint orbit O m = {/i = Coad 9 (yU ) | g G G} = GjG m where G w is the stabilizer 
of /i e 0*. 

2.1 Colored light rays 

Specializing this construction to the case G = SE(3), with C = p 2 and p > 

together with C = 0, we can choose /io = (0, po) and po = (p, 0, 0). 

2 Thc coadjoint representation, Coad, is defined by Coad g /i = /io Ad ff -i where Ad is the adjoint 
representation. 
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The associated 1-form on the group reads then 3 



w = p(u, dx) (2.3) 

where R = (u, v, w) is viewed as an orthonormal, positively oriented, basis of R 3 . 
Its exterior derivative, a = dzu, retains the form 

a(Sg, S'g) = p [<5u, <J'x> - (S'u, 5x>] (2.4) 

for all Sg, S'g G T g SE(3); this 2- form clearly descends to the spherical tangent bundle 
STIR 3 = S 2 x K 3 of R 3 described by the pairs (u, x). 
Computing the kernel of the latter 2-form yields 

(5u,5x) G ker(a) ( * U = ° (2.5) 
v ' y ' I ox = au v ' 

with a G R. We recognize in (|2.5J) the foliation defining the equations of the 
geodesies of Euclidean space (E 3 , ( ■ , ■ )), or, in the context of geometrical optics, 
the light rays in vacuum. 

The first-integrals £ = (q, u) of the foliation (|2.5J) where 

q = x — u(u, x) (2.6) 

parametrize the manifold O M0 = TS 2 of light rays of color p, see (|2.1jl . (This 
invariant, p, of the coadjoint representation of E(3) has been coined "color" in 
as 2irh/p may be interpreted as the Euclidean "wavelength" of the light rays.) The 
symplectic form, u>, of the latter manifold is such that o = (SE(3) — > O^ )*u, viz., 

a(5g, S'g) = lu(S£, S'O = p [{Su, S'q) - (S'u, 5q)] (2.7) 

or, equivalently, 

u = d6 & 6 = -p(q,du). (2.8) 

Let us note that the SE(3)-coadjoint orbit O w is an E(3)-coadjoint orbit, as well. 

3 We write (I2.3|l as a useful shorthand for zu — pSijifdxi , where i,j = 1,2,3; the Einstein 
summation convention is being understood. 
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2.2 The Fermat equations 

In order to describe light rays in a refractive medium of index n G C 1 (IR 3 , M >0 ), let 
us modify the 1-form ()2.3|) according to Fermat's prescription, w -w nw, and start 
with the new 1-form 

tjj = pn(x) (u, dx) (2-9) 

on the bundle STM 3 . 

Again, the characteristic foliation of a = dw should lead to ordinary differential 
equations governing light rays in such a medium. Indeed, (Su, Sx.) G ker(a) iff 
(S(nu), S'x) — (S'(n)u + nS'u, Sx) + A(u, S'u) = for all <5'u, S'x G M 3 , where A is a 
real Lagrange multiplier. Redefining a = X/n, we get 

, r r , . . s f 5(nu) = agradn , . 

(5u,5x)Gker(a)^| ^ = ^ (2.10) 

with a G R, i.e., a 1-dimensional foliation which yields — if we put a = St where t 
is now arc- length — Fermat's equations of geometrical optics in an isotropic medium 
of refractive index n. Notice that the system (|2.1()jl is independent of the color p. 

2.3 The spinning and colored Euclidean coadjoint orbits 
2.3.1 The manifold of circularly polarized light rays 

Apart from the trivial coadjoint orbit and the 2-spheres characterized by the in- 
variants C = and C" = s 2 with s > 0, there exists, most interestingly, another 
class of SE(3)-coadjoint orbits defined by the invariants C = p 2 , with p > (color), 
and C = sp where s ^ stands for spin. The orbit passing through /^o = (^cb Po) 
where £ Q = (s, 0, 0) and p = (p, 0, 0) is again O w = TS 2 and is endowed with the 
symplectic structure coming from the 1-form (|2.2j) on the group SE(3) which reads 

w — p(u, dx) — s(v, dw). (2-11) 
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This 1-form is, by construction, SE(3)-invariant. The associated momentum 
mapping, see J : SE(3) — > O w C sc(3)* : (R, x) i— > (£, p) is actually given by 
^^zfl 1 ) = </(flO ' Z where g i— > is the fundamental vector field associated with 
Z G se(3). li Z = (u, 7), we get <5z(u, v,w,x) = (wxu,wxv,wxw,wxx + 7). 
We readily obtain zu(Szg) = p(u, uj x x + 7) — s(v, a; x w) = (£, u>) + (p, 7), hence 

= xxp + su 

(2.12) 

p = pu 

allowing us to interpret £ as the angular momentum, su as the spin (or polarization) 
vector and p as the linear momentum of the light ray We call helicity the sign of 
the spin invariant, x = sign(s). 

Note, en passant, that the union of two SE(3)-coadjoint orbits defined by the 
invariants (p, s) and (p, —s) is symplectomorphic to a single E(3)-coadjoint orbit. 

An intermediate stage between these classical models and their geometrically 
quantized version is prequantization [31]. Here, the latter construction would restrict 
spin to be half integral, s G ^Zh. To describe spinning light rays, we will naturally 
put s = x^- Euclidean coadjoint orbits with x = +1 (resp. x — ~ 1) describe 
right-handed (resp. left-handed) circularly polarized light rays [3U 122]- For the 
sake of completeness, let us mention that one often writes p = kh where k may be 
interpreted as the (Euclidean) "wave number" of the light rays. 

Straightforward computation yields d(v,dw) = Surf, i.e., the surface element 
of the 2-sphere described by u = v x w, namely Surf (5u, 5'u) = (u, 5u x 5'u). The 
exterior derivative a = dvj takes on the form 

a(5g, 8'g) = p [(5u, <5'x) - (5'u, 5x)} - s(u, 5u x 5'u). (2.13) 

for all Sg,S'g G T 9 SE(3); this 2-form again descends to STIR 3 . The characteristic 
foliation of the latter 2-form is, verbatim, given by ()2.5|) : the spinning light rays in 
vacuum are nothing but the Euclidean geodesies. (As will be shown in the sequel, 
things will change dramatically for such light rays in a refractive medium.) 
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The manifold O^ = TS 2 of spinning light rays is, just as before, parametrized 
by the pairs £ = (q, u) and endowed with the "twisted" symplectic 2-form, u, viz., 

a(5g, 6'g) = co(5^ 6'£) = p [(8u, 8'q) - (5'u, 5q)\ - s(u, 5u x 5'u). (2.14) 

2.3.2 A noncommutative wave plane 

A new phenomenon then appears, which we link to the previously introduced twisted 
symplectic structure, namely noncommutativity of each fiber T U S 2 viewed as the 
"wave plane" u -1 orthogonal to the direction u of the ray. Indeed, let us define linear 
coordinates = (vj,q) with Vj G M. 3 for % = 1,2 in that plane. Straightforward 
computation of the Poisson bracket of q\ and q 2 , with respect to the symplectic 
structure ui in (|2.14j) . yields {91,52} = — uj l (dqi,dq 2 ) = (l/p 2 )(£,v 1 x v 2 ) where 
the angular momentum £ is as in (|2.12j) . If v x and v 2 form an orthonormal basis 
of u^- such that vi x v 2 = u, then 

{9i,? 2 } = 4 ( 2 - 15 ) 

implying noncommutativity of each (Lagrangian) space u -1 consisting of all rays 
parallel to the direction u. 

2.3.3 Euclidean polarized light rays as stationary massless states of the 
Poincare and Galilei groups 

Let us discuss here, for the sake of completeness, the relationship between this purely 
Euclidean model of spinning and colored light rays and the classical models of mass- 
less particles pertaining to either the relativistic or the nonrelativistic framework. 

• The spin-s, massless, coadjoint orbit of E(3, 1) = 0(3, l) xR 4 , the connected 
Poincare group, has the topology of the cotangent bundle of the punctured light cone, 
described by (x, p) G T(M 3 \{0}), and is endowed with the symplectic 2-form ui given, 
in a Lorentz frame, by 

£(5(x, p), 5'(x, p)) = (Sp, 6'q) - (5'p, 5q) - j^(p, 6p x S'p) (2.16) 

HpII 
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corresponding to Equation (17.145) in see also [TBI EH]- The Lie subgroup 
(R, +) of time translations acts in a Hamiltonian way on it. Its infinitesimal action 
is given by 5 e (x, p) = e(— cp/||p||, 0) where eel and c denotes the velocity of light 
in vacuum. The associated moment map is the energy E : (x, p) i— > c||p||, namely 

The submanifold %e '■ STR 3 T(R 3 \{0}) defined by a fixed positive energy 

E = cp (2.17) 

where the constant p > stands for "color", is indeed the presymplectic manifold 
(STIR 3 , a) given by ()2.13|) . where u is the direction of p. This presymplectic structure 
is preserved by the group A(3, 1) = SE(3) xic E(3, 1) . (The "Aristotle group", 
A(3, 1), is the centralizer of time translations in the Poincare group E(3, 1).) The 
Marsden-Weinstein (MW) reduced symplectic manifold T(R 3 \ {0})//(R, +), i.e., 
the stationary relativistic states of energy ([2. 17)1 . is then symplectomorphic to our 
coadjoint orbit (TS 2 ,io) of the Euclidean group SE(3) = A(3, l)o/R, viz., 

i* E w = n*u (2.18) 

where n : (x, p) i— > (q = x — u(u, x), u = p/||p||), and to is as in (|2.14j) . 

• The (massless) coadjoint orbit of Gal(3, l)o = SE(3) x M 4 , the connected 
Galilei group, with Casimir invariants p (color) and s (spin) is plainly symplecto- 
morphic to {TS 2 x TTR, to) described by the quadruples (q, u, t, E), where 

Q = uj — dE A dt (2.19) 

and to is given by (|2.14)1 : cf. Proposition (14.53) of |34j . Time translations, (R, +), 
act in an Hamiltonian way via (q, u, t, E) i— > (q, u, t + e, E) where e G WL. The 
associated moment map is clearly given by the energy (q, u, t, E) \— > E. 

The submanifold je '■ TS 2 xl^ TS 2 x TM + of Galilei massless states of fixed 
energy 

E = const. (2.20) 
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is endowed with the presymplectic 2-form j* E u). Again, its structure is preserved 
by A(3, 1) . (Note that A(3, 1) = E{3, l)nGal(3, 1) C GL(5,R).) The MW-reduced 
symplectic manifold (TS 2 x TR)//(R, +), i.e., the nonrelativistic stationary states 
of energy ()2.20|) . is, again, symplectomorphic to our SE(3)-coadjoint orbit, viz., 



where 7T : (q, u, t) 1— > (q, u), and uj is as in (j2.14j) . 

We have, hence, shown that the connected symplectic manifolds of Euclidean 
spinning light rays, naturally arise as the MW-reduced manifolds of stationary states 
of given energy, E = const., for the massless spinning coadjoint orbits of the con- 
nected Poincare — as well as Galilei — group. 4 

3 A general framework for spinoptics 
3.1 The Fermat prescription revisited 

Let us recall how the Fermat equations describing light propagation in an isotropic 
and inhomogeneous medium can be interpreted as those of the geodesies of a metric 
conformally related to the flat spatial Euclidean metric; see also [TTH ITT] . 

Returning to the expression ()2.9j) which we rewrite as w = pn(x) 2 (n(x) _1 u, dx), 
we introduce the new, curved, metric g = n 2 ( ■ , • ) on M 3 and the g- unitary vector 
U = n~ 1 u; we also put X = x to keep the notation coherent. With these prepa- 
rations, (M = M 3 ,g) becomes a Riemannian 3-manifold while its spherical tangent 
bundle STM = {rj = (U, X) e TM \ g(U, U) — 1} becomes endowed with the 1-form 



which (up to an overall constant factor p) stems from the canonical 1-form of T*M 
via the metric g. 

4 The two-component £?(3)-coadjoint orbits are reduced MW-manifolds of the orthochronous 
Poincare and Galilei groups. 



* — ~ ^* 

J E UJ = 7T U 



(2.21) 



zn=pg(U,dX) 



(3.1) 
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We already know that the characteristic foliation of a = dw yields the Fermat 
equations, see (|2.1()jl . Let us briefly recall why this foliation also provide us with the 
geodesic flow associated with the Fermat metric g = gij{X) dX l ® dX\ where 

gij (X) =n(X) 2 8 ij (3.2) 

with i,j = 1,2,3. (Note that the metric g = ni(X) 2 5ij dX % Cg> dX^ would readily 
enable us to deal with some special anisotropic media JD]-) Call V the Levi-Civita 
connection of (M, g) and its local components. We will denote by <i v the exterior 
covariant derivative of tensor fields, e.g., locally, d^U k = dU k + Y\-dX % U^ . Then, 
the 2- form a = dw = pgij(X) d^U 1 A dX^ writes alternatively 

a(5rj, 8'r,) = p [g(5 v U, 5'X) - g(5' v U, 6X)] (3.3) 

for all Sr), S'rj e T v STM. Its characteristic foliation is integrated by the geodesic 
flow. Actually, 5(U, X) G ker(a) iff g(6 v U, 5'X) - g(5 /V £7, 5X) + Xg(U, 5 ,y U) = for 
all 5 /V £7, 5'X e TM, where A G R. We hence obtain 

5(£/,*)eker(<7)^{ 5 ] U X Z ° xu (3-4) 

with AgI, which we recognize as the geodesic foliation for (M, g). 

3.2 Spinoptics &; minimal coupling to a Fermat metric 

Let us now tackle the geometric description of spinning light rays in an arbitrary 
Riemannian, orientable, 3-manifold (M, g, vol g ). As we will see, this novel approach 
allowing for a complete treatment of the geodesic deviation of spinning light rays in 
a generalized Fermat metric borrows much from general relativity, namely from the 
Papapetrou-Dixon-Souriau equations of motion of test particles in the gravitational 
field, see |S2l ECU ES] and, e.g., [HI ESI d ES! • 

3.2.1 Minimal coupling 

The procedure involved is already known as "minimal coupling" to a classical ex- 
ternal field. In our context, it will simply consist in considering, instead of the 
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Euclidean group (viewed as the trivial SO(3)-principal bundle over E 3 ), the bundle, 
SO(M) — > M, of oriented orthonormal frames of M and in the replacement 

(•,-)~*g & d^d v (3.5) 

into the 1-form (|2.11|) with invariants p and s. With the above notation, this enables 
us to generalize (|3.1|) with the following 1-form 

w = P g{U,dX) - sg(V,d v W) (3.6) 

on the SO(3)-principal bundle SO(M) — > M described by g> = ((£/, V, where 
([/, V, W) is a g-orthonormal basis of T X M such that vol g (Z7, V, W) = 1. 

Taking advantage of the above observation that the characteristic foliation of 
a = dru yields the equations of light rays, we contend that the equations of "geo- 
metrical spinoptics" based on (M, g), for color p and spin s = \h (with \ = ±1), 
are associated with the foliation ker(u) we are now ready to determine explicitly. 

3.2.2 Notation and miscellaneous formulae 

We denote by j(5'X) the g-skew symmetric (cross-product) operator of TxM defined 
by g(6X,j{6'X)6"X) = vo\ g (6X,6'X,6"X). Putting, e.g., V = g(V) = g(V, ■ ), we 
have 

U = j(V, W) <^> ]{U) = WV - VW. (3.7) 
The curvature, R, of the Levi-Civita connection V of (M, g) is defined by 

R(6X, 5'X)5"X = 5 V 5' V 5"X - 5 ,V 5 V 5"X - [5, 8') V 8"X (3.8) 

where X \— > [5, 8']X is the Lie bracket of the vector fields X i— > SX and X i— > S'X. 
Its local expression is given by Rfj k de = ViVjdk — VjVidk, where dk = d/dX k , for 
= !)•••) 3. 

The Ricci tensor Ric(5X, <5'X) = Tr(5"X ^ R(5"X, SX)8'X) has local expres- 
sion R jk = R\ jk . 
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Now, since dim(M) = 3, we have 

R\ 3k = - (R lk 5} - R jk 5f + R\ g ik - R\ g jk ) + l -R (g ik S\ - g jk Sf) (3.9) 

where R = Rijg^ is the scalar curvature. 
If Q = j(U) is as in (|3.7|) . we find 

g(V, R(5X, 5'X)W) = ^g(5X, R(Q)6'X) (3.10) 

where the operator R(Q) is given, via ()3.9|) . by the g-skew symmetric operator 

R(n) = -2(Ricfi + fiRic) + RQ. (3.11) 

The scalar function R(Q, fi) = — Tr(i?($7)$7) will also be needed; from (|3.11j) and 
Q 2 = UU — 1, we get the remarkable expression 

^R(fi, Q) = Ein(£7, U) (3.12) 

where Ein = Ric — \Rg is the Einstein tensor of the metric g. 

For the conformally flat metric ()3.2j) . we readily find the Christoffel symbols 

r* = - (d % n 8* + djn 5* - d e n 5%) , (3. 13) 

the Ricci tensor 

2 11 

Rij = —dindjn didjin An5ij, (3-14) 

Th n n 

where An = S^didjU, and the scalar curvature 

R= _|U n || 2 _ _An. (3.15) 

n 4 " " n 3 v ; 

3.2.3 The general system for spinoptics 

Let us work out the expression of the 2-form a = dm on SO(M) where w is as 
in (EU. 
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We find, remembering ()3.8j) . 

a{8g,8'g) = 8{w{8'g)) - 8'{w{8g)) - w{[8,8']g)) 
= p[g(5 v U,5'X)-g(5' v U,5X)] 

-s [g^V, 5 ,V W) - g{5 ,X7 V, 5 V W) - g(V, R(8X, 8'X)W)] 

for all 8g, 8'g G T g SO(M). Now, using the closure formula UU + VV + WW = 1, we 
readily get g(<5 v V, 8' V W) = g(<5 v V, [1717 + VF + IW]<5 ,V W0 = g(8 w V } UU8' W W) = 
g(8 v U, VW8' V U) since (U, V, W) is a g-orthogonal frame. We then deduce that 
g{8 v V, 5' V W) - g{5' v V, 5 V W) = g{5 v U, [VW - WV}5' V U) = -g(5 v U,j{U)5' v U) = 
vol g (£7, 5 V £7, 8' V U) in view of (jSZZD. At last, we get 

a(8g,8'g) = p [g(5^U, S'X) - g(5' v f/, SX)] 

--sg(8X,R{tt)8'X) (3.16) 
-svol g (f/,(5 v f/,5 /V ?7) 

with the help of ()3.10|) and the shorthand notation Q = j(U). The 2-form ()3.16|) 
turns out to descend, again, to the spherical tangent bundle STM, described by 
7] = (U,X). Still denoting a that 2-form on STM, we compute its kernel. 

Introducing, just as before, a real Lagrange multiplier, A, for the constraint 
g(U, U) = 1, we readily find that 8(U,X) G ker(a) iff p5 v U + \ sR(tt)5X = and 
p8X + s]{U)5 v U = XU. This entails p8X - s 2 / {2p)VlR{tt)5X = XU, prompting the 
Ansatz 8X = all + (3s 2 QR(Q)U, for some a, (3 G IR still to be determined. We thus 
get p 2 (aU+f3VLR{ti)U)-\ s 2 aVtR{n)U-\ s 4 /3fi \R{tt)VLR{tt)} U = pXU. Recall that, 
if A, B are g-skew symmetric operators of TxM, then ABA = ^Tr(AB)A. This en- 
ables us to compute the above bracketed term, viz., R(Q)QR(Q) = — ^R(Q, Q)R(Q), 
and to find a = X/p and (3 = X/{2p{p 2 + | s 2 R(Q, Q)). Invoking (|3.12|) . we end up 
with the spinoptics system governing the trajectories of spinning light rays in a 
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Riemannian 3-manifold (M, g), viz., 




pS v U 



lsR{Q)5X 



5(U,X) e ker(cx) 



s 2 QR(n)U 



(3.17) 



5X 



a U + 



2 [p 2 + s 2 Ein([7, U)} 



with Q = j(U), and aeR. 

In the (3 + l)-dimensional setting of general relativity, a similar system (the 
Papapetrou-Dixon-Souriau equations) would describe geodesic deviation and spin 
precession of spinning test particles in a gravitational background field. 

We stress that, in spinoptics (as well as in general relativity [SH]), "velocity" 5X 
and "momentum" U fail to be parallel (see (|3.17jl ). We will see how this phenomenon 
gives rise to subtle physical effects such as the Optical Hall Effect. 

3.2.4 Fermat spinoptics 

As an illustration of the preceding results, let us write the exact spinoptics equations 
specialized to the Fermat metric (|3.2|) associated with a refractive index n. Just as 
before, we denote by p(= kh) the color and s(= xty the spin of the model. Although 
this could be deduced from the very general system (|3.17j) . we choose to simply start 
from the 1-form ()3.6|) as this procedure actually yields all parameters adapted to 
the model in a straightforward fashion. 

We begin here with (M, g) = (IR 3 , n 2 ( ■ , ■ )) where n G C 2 (1R 3 , M >0 ), and, upon 
defining the Euclidean frame (u, v,w) = (nlI,nV,nW), express the 1-form ()3.6j) 
as a new 1-form on SE(3) which is parametrized by g = ((u, v,w),x) as in Sec- 
tion^ Using the expression (J3.13|) of the Christoffel symbols, we readily notice that 
g(V, d v W) = (v, dw) + ra(x) _1 [(v, dx) dn(w) — (w, dx) dn(v)]. Introducing now the 
"velocity" 



1 



(3.18) 



v = — 
n 



and its gradient 



g = grad(w), 



(3.19) 
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we find the new spin term g(V, d v W) = (v, dw)—n [(v, dx) (w, g) — (w, dx) (v, g)] = 
(v, dw) +n((vxw)x g, dx). 

This entails, that w = pg(U, dX) — sg(V, d?W) retains the following form 

w = (p,dx) - s(v,dw) (3.20) 

where 

p = n(x)(pu + sg x u) (3-21) 

can be consistently interpreted as the spin-dependent "momentum" of the system. 

Easy computation gives the exterior derivative a = dw of the 1-form ()3.20|) . 
namely 

a(6g, 8'g) = (8p, 5'x) - (6'p, 8x) - s(u, 5u x 5'u). (3.22) 

for all 5g, 5'g G T g SE(3), see ()2.13|) . This 2-form descends to STR 3 and has, gener- 
ically, rank 4; computing its kernel needs some more effort. 
Let us denote by 

Vg = §J (3.23) 
ax 

the (symmetric) second derivative of the velocity v. 

We find, using (|3.21j) . and in the same way as before, (5u, 5x) e ker(cr) iff 
5p + ng(p, Sx) — ns(Vg)u x 5x = and pnSx — nsg x 5x + su x 5u = Au where 
A G M. is a Lagrange multiplier. Taking the cross-product of the latter equation 
by u yields s5u = u x n(p5x — sg x 5x), which can be inserted into the former 
equation with the help of (|3.21|) . In doing so, using the following partial result 
5u + (s/p)g x 5u = [(np/s)u + ng x u + (sn/p)g(g, u)] x 5x, we end up with 



u x 5x + n 2 s||g|| 2 u x 5x - sn j(u) Vg + Vg j(u) Sx = (3.24) 



s 

where, see ()3.7j) . j(u) : 5x i— > u x 5x is the Euclidean cross-product operator. We 
still need to compute the last anticommutator in the above equation; it is given by 
the general formula j(u) H + H j(u) = j(—Hu + Tt(H)u) for any u G 1R 3 and any 
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symmetric H £ L(M 3 ). Equation ()3.24j) reduces then to 



5x = p 



r p 2 n 2 



u + n 2 s||g|| 2 u + sn(Vg)u — sn Tr(Vg)u 



for some (3 E M. 

We can finally write the system defining the kernel of our 2-form as 

(5u, £x) e ker(cr) 

t 

P 



s5u 
5x 



U X 1 j 

V p 

-,2 



5x 



a: 



with aei and where 



a = 1 + 



2 „ || 2 -^div(| 



p* p 
We have, equivalent ly 5 



au + — V u g 



,, , ll2 s 2 ■ 
IH — 5-||gradf|| -uAi; 



P 



p- 



5p 



5x 



(<5p, 5x) e ker(cr) 
s dp. 



-(p,5x g + 
f pox 



s 

1 + -J( 

p 



-1 



x 5x 



f s 2 dg 
pr ox 



s 
P 



(3.25) 



(3.26) 



(3.27) 



(3.28) 



with a G M, in terms of the natural variables, viz., momentum p and position x. 

This rather complicated system defining the equations for the trajectories of 
spinning light rays with color p and spin s in a refractive medium of index n = 1/v 
constitutes the novel differential equations of Fermat spinoptics, up to reparametriza- 
tion. 

We clearly recover from (|3.28|) the original Fermat equations (I2.10J) in the spin- 
less case, s = 0. 



5 Note that (1 + j(z)) 1 = (1 + ||z|| 2 ) -1 (1 - j(z) + zz), where z = (z, • ), for all z S 
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Let us finish this Section by highlighting the relationship of our foliation for 
Fermat spinoptics to recent work of Onoda, Murakami and Nagaosa [S3] . Neglecting, 
in our system (|3.28|h all terms involving second derivatives Vg of the velocity v, and 
all quadratic terms ||g|| 2 , we end up (choosing a parameter defined by a = 1), with 
the system 

Sp — <p,<Jx)g 

<5x = p gxp 

V 

which, up to notation and reparametrization, exactly matches the first two Equa- 
tions Of Motion put forward in [S3] . The EOM (J3.29j) provide, hence, a linearization 
of our system ()3.28|) around g = 0. 

4 Spin Snell-Descartes' laws & Optical Hall Effect 

We will, as a first test of our approach, establish the spinoptics version of the Snell- 
Descartes laws generalizing those of plain geometrical optics. 

Consider the simplest case of a planar interface separating space into two regions 
(Mi, Hi), resp. (M 2 ,n 2 ), where Mi = {x G M 3 |(n, x) < 0} has refractive index 
Hi = const., resp. M 2 = {x 6 I 3 ( n ) x ) > 0} has refractive index n 2 = const., 
where n is a unit vector, orthogonal to the interface (and pointing toward M 2 ), 
characterizing the optical device. 

Wishing to describe the laws of reflection and refraction, in geometrical terms, 
namely the scattering of spinning light rays by this device, we will resort to a theory 
developed by Souriau, namely "symplectic scattering" [33]. See also [2SJ for some 
further developments. 

4.1 Symplectic scattering 

Symplectic scattering should be thought of as the classical counterpart of unitary 
scattering of quantum mechanics or quantum field theory. Classically, what is pre- 
served by a scattering diffeomorphism is the basic structure of the theory, namely 
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the symplectic structure of the manifold of classical states, whereas quantum me- 
chanically, the scattering S"-matrix has the property of preserving the fundamental 
structure of the space of quantum state vectors, namely the Hilbertian structure. 

Given symplectic manifolds (A4\,uji) of "in" states, and (Ai 2 ,u 2 ) of "out" 
states, we assume that a scattering process is given by a local symplectomorphism, 
viz., a local diffeomorphism 

S : Mi -> M 2 such that k>i = S*u 2 . (4.1) 

Such mappings being far from unique, we have to take into account the geometric 
features of the scattering device to try and find a unique symplectomorphism, S. 

In most cases, the "in" and "out" manifolds are Hamiltonian G-spaces 6 (e.g., 
coadjoint orbits) of some Lie group G, for instance a group of space(-time) auto- 
morphisms; they represent the free asymptotic states of the system. The scattering 
device reduces the original symmetry to a Lie subgroup H C G whose action is 
assumed to intertwine the symplectomorphism S, that is 

S o h Ml = h M2 o S (4.2) 

for all h G H. If Z G f), where f) is the Lie algebra of H, we readily find from (|4.2jl that 
the associated fundamental vector fields are S'-related, Z M2 = S*Z Ml . This entails, 
via flUJ) and flOD, that u x {Z Ml ) = S*(u 2 )(Z Ml ) = S*(u 2 (S*Z Ml )) = -S*{d{J 2 -Z)) 
where J 2 is the moment map of (M. 2 ,u 2 ,G). At last u^Zj^n) = —d(S*(J 2 ) ■ Z) for 
all Z G f) which, if Ai± and A4 2 are connected, enables us to write the conservation 
law 

Ji|& = S%7 2 |f)) (4.3) 

that plays a central role in the determination of the sought scattering mapping, S. 

6 These are symplectic manifolds (M, ui) equipped with a G-action g i— ► for which g* M io = u> 
and a momentum mapping 34 , i.e., a globally defined mapping J : M. — > g* such that there holds 
v(Z M ) = ~d(J ■ Z) for all Z e g. 
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4.2 Scattering of spinning light rays 

Let us now turn to the effective computation of the scattering mapping of light rays 
with color p and spin s = ±h by the previously introduced interface separating two 
media of constant, unequal, refractive indices n\ and n 2 - 

From now on {M.\,UJ\) and (M. 2 ,uj 2 ) will represent SE(3)-coadjoint orbits re- 
spectively characterized by the invariants 

Ci=pl C[=p lSl , & C 2 =p 2 2 , C' 2 =p 2 s 2l (4.4) 

where p\ = pn\ & p 2 = pn 2 with p > and s±, s 2 E {+h, —H}. We handle, in this 
manner, all helicities at the same time. 

The canonical 2- forms on M. i = TS 2 are given by ()2.14j) and read now 

Wi(<*&, = Pi [(^ u i ; - (S'ui, 5qi)) - Si(ui, 5ui x 6'vu), ( 4 -5) 

for i = 1, 2. 

Incoming light rays, i.e., hitting the interface in M±, constitute a submanifold 
of whereas light rays refracted in M 2 form a submanifold of Ai 2 . Moreover, 
reflection will be dealt with by considering A4 2 = M. manifold, whose sym- 

plectic 2-form u 2 is defined by p 2 = p\ (since light bounces back in half-space M\ 
with index %) and s 2 . 

As to the symmetry group of the optical interface, it is clearly given by the Lie 
subgroup 

H = {(A, c) G SE(3) I An = n, (n, c) = 0}, (4.6) 

hence H = SE(2) C SE(3). 

We are now ready to implement (|4.1jl . ()4.2)1 and (|4.3)1 . 
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4.2.1 Conservation laws 

The if -momentum mapping of (TS 2 ,uj) is the restriction J|f) of the Euclidean mo- 
mentum mapping, J. We find [31] that J(q, u)|f) = (L, P) is of the form 

L = (n,£) 

(4.7) 

P = n x p 

where J(q, u) = (£, p) is as in (J2.12jh 

If we put (q 2 ,u 2 ) = ^q^iix), the conservation law (|4.3|) reads 

(n, qx x pi + s x ui) = (n, q 2 x p 2 + s 2 u 2 ) (4.8) 
nxp! = n x p 2 (4.9) 



where 



Pi=PiUi (4.10) 



for i = 1, 2. 

Equation ()4.9|) readily implies 

p 2 = Pi + An (4.11) 

where A is some a smooth function of (qi,ui); taking into account the Euclidean 
invariants C\ = ||pi|| 2 — p\ and C 2 = ||p2|| 2 — v\t see f!4.4|) . already insures that A 
depends on 

a=(n,pi) (4.12) 

only, via 

A 2 + 2aA + d-C 2 = 0. (4.13) 

Note that, with our orientation, incoming rays are such that a > 0. We readily find 
the explicit expression 7 



-a + y/a 2 + C 2 — Ci (refraction, rii ^ n 2 ) 

A = <j (4.14) 
—2a (reflection, ri\ = n 2 ) 

which will be used in the sequel. 

7 If Ci > C2, then a 2 + C2 — C% > must furthermore hold true; if the latter condition is not 
satisfied, total reflection occurs. 
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4.2.2 The scattering symplectomorphism 

Taking advantage of (|4.11|) . we now seek the diffeomorphism S : (qi, Ui) i— »■ (q2, u 2 ) 
starting from the general Ansatz 



(4.15) 



(4.16) 



q2 = qi + A*Pi + z/n + x pi 
P2 = Pi + An 

where A is given by (|4.1H|) and /i, is, g are otherwise arbitrary functions of (qi, Ui). 

• From (|4.15jl . (|4.8jl . (|4.4j) . together with C'/C = s/p, we immediately obtain 
a(C' 2 /C 2 — C[/Ci) — g(Ci — a 2 ) + XC' 2 /C 2 = 0, or, if n x p x ^ (in the generic case 
of non normal incidence), 

1 

||n x p x | 

In the case of normal incident rays, ||n x pi|| 2 = C\ — a 2 = 0, we must have 
(A + a)C' 2 /C 2 - aC'JCx = with a = y/C[ and A + a = y/C~ 2 (resp. X + a = -y/C[) 
for refraction (resp. reflection). We therefore find 

s 2 = s\ (refraction) 
s 2 = — s\ (reflection) 

which constitute non trivial conditions on the scattering symplectomorphism, S. 

• Take now into account the constraints 



C 2 C[\ C' 2 



(4.17) 



I P. 



Ci & (p h qj) = 0, 



(4.18) 



for all i — 1, 2, to further determine the yet unknown function v. If we put 



z = (n,qi) 



(4.19) 



for the n-component of "position" qi, then (pi,qi) = (P2,q2) = 0, together 
with (J4.15j) imply 



v 



-1 

a + A 



(Xz + n(Cx + aX)) . 



(4.20) 



24 



• Let us use the previous Ansatz (|4.15|) to express that S : £1 i— > £2 is a 
symplectomorphism, namely = 1^1(^2,^2) for all 6£i,5'£i e T^A^, or 



C" 

(5pi,5'qi) - (5'pi,5qi) - ^|(pi,5pi x 5'pi) 



(4.2i; 



(<5p 2 ,5'q 2 ) - (5' P2 ,5q 2 ) - t|(p2,5p2 x 8'p 2 ) 

^2 



for all tangent vectors compatible with the constraints (|4.18j) . 
A tedious calculation shows us that (|4.21j) yields 

= +5X5' z - 5'X5z + a(6\6'fi - 5'X5^) + 5a5'v - 5'a5v + 5X5' v - 5'X5v 

+(Spi, n x Pl )5'(g - C' 2 /C 2 2 X) - (5' Pl , n x Pl )5(g - C' 2 /C 2 2 X) (4.22) 
-((2g + XC' 2 /C 2 2 )n,5 Pl x 5' Pl ) + (Cj/Cj - C' 2 /C 2 2 )( Pl , 5 Pl x 5' Pl ). 

In order to tackle (|4.22j) . we find it useful to introduce spherical coordinates (9,ip) 
on the 2-sphere described by ui = P i/pi = (cos <p sin 9, sin ip sin 9, cos9). 8 
Rewrite as 

= +dX Adz + adX A dp, + da A dv + dX A dv 

+d sin 2 A dg - C 2 /C 2 tZA) (4.23) 
-{2q + XC'JCDd cos sin 9 d9 A dtp + {C[/C 2 - C 2 /C 2 )Cf /2 sin 9 d9 A dtp. 

From (jUinD we get dX = -Xda/(a + X), while (P~T2|l yields rfa = -y/Chsm 9d9. 
We then obtain the partial expression dX A dz + ct<iA A dp + da A dv + dX A dv = 
—l/(a + X)da A d(Xz + aXp — av) = —l/(a+ X) 2 da A d((C 2 — C\)z + C 2 a\i) with 
the help of ()4.20|) . Some more effort is needed to finally transcribe ()4.23j) as 







da A 



d 



[C 2 - C x )z + C 2 ap) 
(a + X) 2 

r fC 2 C[ 
°i I ^2 ~~ 7^2 



\ rv sii rv 
A ^2 °2 1 u l 



+ ) 



c 2 c 2 



(a + X)C 2 C 2 d 



C 2 ^ (d-a 2 

aA cf + A ^^TA) 



(4.24) 



8 We have (ui, Sui x J'ui) = sin 
and also (n, Ui x dui) = sin 2 9 dip. 



'tp-S'eStp), and (n, <Jui x 5'm) = costf sin6(60S'cp-S'eScp), 
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This readily implies 



C 2 - d + C 2 J^ = (4.25) 
oz 



and, with the help of ()4.13|) . also gives 



which leaves us with 



C 2 a 



where jtt = /(a) is an arbitrary function of a. 

• So far, all four functions \,{i,v,q have been determined by (|4.14j) . ()4.27|) . 
(|4.2()jl and (|4.16J) . up to an arbitrary function ju. Let us show that, indeed, p, = 0. 

Returning to the expression (j4.15|) giving the scattering mapping which we 
write, for convenience, S : (qi,Pi) (q 2 = qi+^iPi+^in+^inxp!, p 2 = Pi+Ain), 
its inverse S' 1 : (q 2 , p 2 ) h-> (qi = q 2 + /i 2 p 2 + z/ 2 n + £ 2 n x p 2 , p x = p 2 + A 2 n) is 
such that A 2 = — A x (where a 2 = «i + Ai), /i 2 = — /ii, ^ 2 = Ai/ii — z^i and £> 2 = — £>i. 
These relationships implement the principle of ray reversibility. 

From the definition (j4.19|) . we get z 2 = Z\ + a.\\i\ + z/i and find, resorting 
to (jQ7|l . that /ii + /i 2 = [Cia|/(a 2 ) + C 2 a?/(ai)] /Cia| = iff / = 0. 

We obtain, at last, 

,= ^f. (4.28) 
C 2 a 

We have thus completed the explicit determination of the scattering symplecto- 
morphism by the plane interface separating two refracting media of constant indices. 9 
Let us collect and present the above findings in a new guise where the scattering 



9 The ff-equivariance (|4.2[> of the unique symplectomorphism (|4.15|) . with (|4.14|) . I|4.28|) . I|4.2U[) 

and (JUSJl, can be directly checked to hold, the H = SE(2)-action on M. = TS 2 being given by 
}im(<1, u ) = (Aq + c - Au{Au,c),Au), where h = (A,c) G H (see (14.6^ . 
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mapping, S, is uniquely given by 1)4.15)1 with 1)4.17)1 and 

A = 



(n, Pi) + y/C 2 - ||n x pi| 

2(n,pi) 

(n,qi) 



(refraction) 
(reflection) 



C 2 



(n, Pi) 



(4.29) 



Ci (n,qi) 
C 2 (n,pi) 
1 

||n x pi| 



-(c> 




\c 2 




; as in 


(ED- 



(n,Pi) 



a 



•A 



4.3 Snell-Descartes' laws of spinoptics &; Optical Hall Effect 

Introducing the angle of incidence 9\ (resp. the scattering angle 9 2 ) between pi 
(resp. P2) and n, we easily infer from (J4.9)) the law of refraction obeyed by the 
direction of light rays, namely pisin^i = p 2 sin6 2 . As for the law of reflection, 
(14.14)1 already yields the mirror transformation pi 1— > p 2 = pi — 2n(n, pi). 

4.3.1 The generalized Snell-Descartes laws 

Summing up, and taking into account the specific result ()4.17j) . we write the Snell- 
Descartes laws of spinoptics as 

n 2 sin #2 = n\ sin^, s 2 = Si, (refraction) 

9 2 — 7T — 9\, s 2 = — Si, (reflection) 

Let us emphasize that these laws must be supplemented — as shown below — by a 
new law which unveils a phenomenon pertaining to geometrical spinoptics, namely a 
transverse shift of the scattered spinning light rays off the plane of incidence spanned 
by n and pi in generic position. 

Choose now as origin, O, of Euclidean space, the intersection of the incoming 
light ray and the interface separating the refractive media, so that qi = 0. From (J4.15)) 

and ()4.29)1 . we get \i = and v = 0. We then obtain q 2 = £>n x p x where 

1 



(4.30) 



Q 



n x p! I 



s 2 , . Si 

— n,p 2 ) (n,pi) 

P 2 Pi 



(4.31) 
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Note that there is no transverse shift in the case of normal incidence. We finally 
obtain the following expression for this transverse shift 

n x p! [s 2 cos 6 2 - Si cos 0i] 

c b~ c Li = ir-^ — ii 432 

ll n x Pill p ?7,i | smt/il 
which is clearly valid for either cases of refraction or reflection. This formula does 
formally agree with an analogous expression proposed by Onoda, Murakami and 
Nagaosa [2] who used quite a different viewpoint. 

Note, however, that the transverse shift (|4.33|) for reflected rays vanishes in 
our framework; see also [SI]. On the other hand, the nontrivial transverse shift for 
refracted spinning light rays, theoretically explained by (|4.33|) in the present context 
of geometrical spinoptics, is indeed a novel phenomenon known as the "Optical Hall 
Effect" and of a great importance in the new trends of experimental optics. 

4.3.2 The special case of left-handed media 

So far, we have been dealing with ordinary dielectric media. Quite interestingly, 
artificial materials enjoying a negative refractive index (with both negative dielectric 
susceptibility and magnetic permeability) have been foreseen by Veselago a few 
decades ago jSH]- These brand new "left-handed" media (LHM) or metamaterials 
are nowadays manufactured in the laboratories and their strange optical properties 
systematically studied. See, e.g., [§]. 

Let us emphasize that our general formalism applies just as well in the presence 
of these LHM. For example, the Snell-Descartes laws, spelled out in Section 14.3.11 
still hold true if n 2 < 0, say. In this case, Equations ()4.30|) account for "nega- 
tive reflection" , while (J4.1U)) shows that the linear momentum p 2 = pn 2 U2 and the 
direction u 2 of a refracted ray are antiparallel. 

One of the striking feature of these metamaterial is that the transverse shift of 
reflected and refracted spinning light rays may vanish identically, namely 

qi - q 2 = iff n 1 +n 2 = 0. (4.33) 
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This characteristic property of metamaterials might provide a purely mechanical 
interpretation of the "perfectness" of LHM lenses that allow to break the Rayleigh 
limit of optical devices [33] . 

5 Conclusion and outlook 

The basics of geometrical spinoptics have been laid to extend, from first principles, 
geometrical optics to spinning light rays. The point of view we have espoused made 
crucial use of Euclidean geometry. By generalizing the Fermat prescription to the 
presymplectic manifolds upstream of generic coadjoint orbits of the Euclidean group, 
we have derived a 1-dimensional foliation governing the trajectories of spinning light 
rays in arbitrary dielectric media. A refinement of the classic Snell-Descartes laws 
readily followed, together with the expression of the local scattering symplectomor- 
phism undergone by spinning refracted and reflected light rays. This enabled us 
to derive a formula for the associated transverse shift, specific to the Optical Hall 
effect. 

A number of queries, triggered by the present study, remain the subject of future 
work; let us mention here some few examples. 

Revisiting the theory of caustics within this new framework would certainly be 
a worthwhile task, in view of the refinement the Snell-Descartes laws governed by 
noncommutativity of the wave plane. 

Quantizing geometrical spinoptics is also serious endeavor. One might profit 
by the fact that prequantization of the symplectic manifold of photonic states, with 
s = ±h, is given by the contact structure on the quotient SO(M) / (ker(w) r\ker(dw)) 
defined by (|3.6p . Another route to quantization might, alternatively, be offered 
by the procedure of conformally equivariant quantization 20 . It would also be 
interesting to see how close to Maxwell theory should such a quantization lead us. 

Another challenging project would be to tackle all polarization states at a single 
stroke by considering that circular polarization states given by the coadjoint orbits 
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of the Euclidean group are, in fact, the building blocks or elementary systems of a 
more elaborate, quantum, theory of spinoptics. 

At last, it seems reasonable to envisage extending geometric spinoptics to the 
case of Faraday-active media by coupling, from the outset, the photon spin and 
the external magnetic field via the color, p, much in the same way as the (anoma- 
lous) magnetic interaction term is introduced via the mass in a general relativistic 
framework [H E3 EH EEE! ■ 
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